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Abstract. We generalize Grover's unstructured quantum search algo- 
rithm to enable it to use an arbitrary starting superposition and an 
arbitrary unitary matrix simultaneously. We derive an exact formula for 
the probability of the generalized Grover's algorithm succeeding after n 
iterations. We show that the fully generalized formula reduces to the spe- 
cial cases considered by previous authors. We then use the generalized 
formula to determine the optimal strategy for using the unstructured 
quantum search algorithm. On average the optimal strategy is about 
12% better than the naive use of Grover's algorithm. The speedup ob- 
tained is not dramatic but it illustrates that a hybrid use of quantum 
computing and classical computing techniques can yield a performance 
that is better than either alone. We extend the analysis to the case of a 
society of k quantum searches acting in parallel. We derive an analytic 
formula that connects the degree of parallelism with the optimal strat- 
egy for fc-parallel quantum search. We then derive the formula for the 
expected speed of fc-parallel quantum search. 



1 Introduction 



The field of quantum computing has undergone a rapid growth over the past 
few years. Simple quantum computations have already been performed using 
nuclear magnetic resonance p6| , ^0[pT| ,P,p|,pT| and nonlinear optics technologies 
iFlP"^- Recently, proposals for more specialized devices that rely on quantum 
computing have also been made p3| . Such devices are far from being general- 
purpose computers, nevertheless, they constitute significant milestones along the 
road to practical quantum computing. 
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In tandem with these hardware developments, there has been a parallel devel- 
opment of new quantum algorithms. Several important quantum algorithms are 
now known f p2p^ , [l7| ,p|.p^||^| . Of particular importance is the quantum algorithm 
for performing unstructured quantum search discovered by Lov Grover in 1996 
[ p"7| . Grover's algorithm is able to find a marked item in a virtual "database" 
containing N items in 0(y/~N) computational steps. In contrast, the best classical 
algorithm requires 0(N/2) steps on average, and O(N) steps in the worst case. 
Thus Grover's algorithm exhibits a polynomial speedup over the best classical 
counterpart. 

Although the Grover algorithm exhibits only a polynomial speedup, it ap- 
pears to be much more versatile than the other quantum algorithms. Indeed, 
Grover has shown how his algorithm can be used to speed up almost any other 
quantum algorithm |is|] . More surprisingly, even search problems that contain 
"structure" in the form of correlations between the items searched over, often 
reduce to an exhaustive search amongst a reduced set of possibilities. Recently, it 
was shown how Grover's algorithm can be nested to exploit such problem struc- 
ture ||. This is significant because NP-hard problems, which are amongst the 
most challenging computational problems that arise in practice, possess exactly 
this kind of problem structure. 

In order to appreciate the full versatility of Grover's algorithm it is important 
to examine all the ways in which it might be generalized. For example, whereas 
the original Grover algorithm was started from an equally weighted superposi- 
tion of eigenstates representing all the indices of the items in the database, a 
natural generalization would be to consider how it performs when started from 
an arbitrary initial superposition instead. This refinement is important, because 
if Grover's algorithm is used within some larger quantum computation, it is likely 
to have to work on a arbitrary starting superposition rather than a specific start- 
ing eigenstate. Similarly, the original Grover algorithm uses a particular unitary 
operator, the Walsh-Hadamard operator, as the basis for a sequence of unitary 
operations that systematically amplifies the amplitude in the target state at the 
expense of the amplitude in the non-target states. However, it is now known that 
this is not the best choice if there is partial information as to the likely location 
of the target item in the database. In such a situation a different unitary operator 
is desirable. Hence, it is important to understand how Grover's algorithm per- 
forms when using an arbitrary unitary operator instead of the Walsh-Hadamard 
operator. 

Each of these refinements have been analyzed in detail independently. Bi- 
ham et al. have considered the case of an arbitrary starting superposition || 
and Grover considered the case of an arbitrary unitary operator fl8|| . In this pa- 
per, we present the analysis of the fully generalized Grover algorithm in which 
we incorporate both of these effects simultaneously. Our goal is to determine 
the exact analytic formula for the probability of the fully generalized Grover 
algorithm succeeding after n iterations when there are r targets amongst N 
candidates. Having obtained this formula, we will recover the Biham et al. and 
Grover results as special cases. We will then show that the optimal strategy, on 
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average, for using the fully generalized Grover algorithm consists of measuring 
the memory register after about 12% fewer iterations than are needed to obtain 
the maximum probability of success. This result confirms a more restricted case 
reported in || . Finally, we show how to boost the success probability and reduce 
the required coherence time by running a society of k quantum searches inde- 
pendently in parallel. In particular, we derive an explicit formula connecting the 
degree of parallelism, i.e., k, to the optimal number of iterations (for each agent 
in the society) that minimizes the expected search cost overall. We then derive 
the expected cost of optimal /c-parallel quantum search. 



2 Grover's Algorithm 

The problem we have to solve is the following. Given a function f(xi) on a set 
X of input states such that 

, , n J 1 if Xi is a target element , . 

}{Xl) ~ \0 otherwise " [ > 

How do we find a target element by using the least number of calls to the function 
f(xi)7 In general, there might be r target elements, in which case any one will 
suffice as the answer. 

To solve the problem using Grover's algorithm we first form a Hilbert space 
with an orthonormal basis element for each input Xi € X. Without loss of 
generality, we will write the target states as \ti) and the non-target states as 
In this paper we refer to the basis of input eigenstates as the measurement 
basis. Let N = \X\ be the cardinality of X. The function call is to be implemented 
by a unitary operator that acts as follows: 

\^)\y)^\x i )\y®f(x i )) (2) 
where \y) is either |0) or |1). By acting on 

f N-r r \ 

£>|ii) +X>l*i>) (|0> |1>) (3) 
with this operator we construct the state 

'N— r r \ 1 

5>io-5>i**> -/= ao>-ii» (4) 
»=i j=i j » 

where the r measurement basis states are the target states and the N — r 
measurement basis states \U) are the non-target states. If we now disregard the 
state -j= (|0) — |1)) then all we have done is to invert the phase of the target 
states. Hence, the operator we have achieved is equivalent to the operator 
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l- 2^; !*,)(*,!. (5) 



although we emphasize that it is not necessary to know what the target states 
are a priori. 

Next we construct the operator Q defined as 

Q = {2\a){a\-l)(l-2^\t i ){t i \\ (6) 

where |a) can be thought of as the averaging state. Different choices of \a) give 
rise to different unitary operators for performing amplitude amplification. In the 
original Grover algorithm, the state \a) was chosen to be 

(7) 




and was obtained by applying the Walsh-Hadamard operator, U, to a starting 
eigenstate \s), i.e., \a) — U\s). Hence, the operation 2|a)(a| — 1, which Grover 
referred to as "inversion about the average", is equivalent to —UI S U> with U 
being the Walsh-Hadamard operator and I s being 1 — 2|s)(s|. 

By knowing more about the structure of the problem we can choose other 
vectors \a) that will allow us to find a target state faster. Techniques for doing 
this are given in M] . 

Fortunately, in order to determine what action the operator Q performs, it 
is sufficient to focus on a two-dimensional subspace. The basis vectors of this 
subspace can be written as 

I*) =$£Li<*il«>l*i>> « 2 = EUIWI 2 ^ 

K) = 7=f (\a)-v\t)) W 

Note that \t) is the normalized projection of \a) onto the space of target states 
and | a') is the normalized projection of \a) onto the space orthogonal to |t). This 
choice of basis makes the calculation easiest. The rest of the Hilbert space (i.e. 
the space orthogonal to \t) and \a!) ) can be broken up into the space of target 
states (St) and non-target states (Si,). We can now write Q as 



Q = cos4>(\t)(t\ + |a')(a'|)+sin0(|t)(a'| - \a')(t\)+I T -I L , = arccos [l - 2v 2 ] 

(9) 

where It and 1^ are the identity operators on (St) and (Sl) respectively. From 
this we can see that Q is just a simple rotation matrix on \a') and \t) and acts 
trivially on the rest of the space. 

An arbitrary starting superposition |s) for the algorithm can be written as 

\ S )=a\t)+[3e lb \a') + \ ( f H ) + \ ( f> l } (10) 
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where the states \4>t) and \<j>i) (which must have a norm less than one if the state 
\s) is to be properly normalized overall) are the components of \s) in (St) and 
(Sl) respectively. Also, a, (3 and b are positive real numbers. After n applications 
of Q on an arbitrary starting superposition \s) we have 



Q n \s) = (acos(n<j)) + f3e lb sm(n</))) \t) + ((3e lb cos(n0) - asin(n<£)) \a') + \<p t ) + (-l) n |<&). 

(11) 

If we measure this state our probability of success (i.e., measuring a target state) 
will be given by two terms. The first term is the magnitude squared of Q n \s) 
in the space St- This magnitude is (<j>t\4>t) and is unchanged by Q. The second 
term is the magnitude squared of the coefficient of \t) which is given by 



g(n)^\(t\Q-\ S )\ 2 

= | a cos(n</>) + (ie lb sin(n</>) | 
= + cos(2n</>) + a/3cos(6) sin(2n<£) 

= + \ \a 2 + (3 2 e 2lb \ cos(2n0 + V) 



(12) 
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where ip = arccos r^r+p^zwi ■ This is the term that is affected by Q and is the 
term we wish to maximize. The total probability of success after n iterations of 
Q acting on \s) is 

p(n,r,N) = (<j> t \<f> t )+g(n). (13) 

Assuming that n is continuous (an assumption that we will justify shortly) the 
maxima of g(n), and hence the maxima of the probability of success of Grover's 
algorithm, are given by the following. 

nj = ^ (-V + 27rj) j = 0,1,2... (14) 
The value of g(n) at these maxima is given by 

9(n j ) = ^±^+\\c?+pe™\. (15) 

In practice, the optimal n must be an integer and typically the n/s are not 
integers. However, since g(n) can be written as 

g(n ±S) = g(n ) - (f> 2 \a 2 + (3 2 e 2lb \S 2 + 0(5 A ) (16) 

around rij and most interesting problems will have v <§C 1 and hence <fi ~ 2v -c 1, 
simply rounding rij to the nearest integer will not significantly change the final 
probability of success. So, we have 

p(n max , r, N) = <?L±$L + I | a 2 + 02 e 2ib\ + { ^ t) _ 0(w 2 } 

as the probability of measuring a target state after n max — rij applications of 
Q. 
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3 Recovering the Special Cases 

As a check on our fully generalized formula for the probability of success after 
n iterations, we attempt to recover the corresponding formulae obtained in the 
analyses of Biham ct al (for a fixed unitary operator and an arbitrary starting 
superposition) Q and Grover (for an arbitrary unitary operator and a fixed 
starting superposition) Q. 

In the case of Biham et al., the starting state is arbitrary but the averaging 
state | a) is given by 



77* E l*> (18) 



N , ,. 



In this case 



I*} =^ET=il*i> (19) 

In the analysis of [|| they use fc(0) and 7(0) to represent the average amplitudes, 
in |s), of the target and non-target states respectively, and Uk and o\ to represent 
the standard deviations of those amplitudes. With some algebra one can see that 
the following relationships connect our notation to theirs: 



a — > k(0)y/f 

— r 

(20) 



f3e lb — ► l(0)VN^ 
(4>t\<t>t) — ► raf 
{4>i\4>i) ► (N — r)ofc 
4> — > LO 
ijj — > 2Re[</>] - 7T 
n — ► t 
n — * T. 

If you substitute these relationships into equations |l2|, [u], and [Tt] you will re- 
produce the results of || . 

The second special case, in which |o) (the averaging state) is an unknown 
normalized vector while \s) is given by 

l s > = I") (ou 

= VT - vS\a')+v\t). [ ' 



was considered by Grover. Hence, a = v, (3 = vl — v 2 and 6 = 0. Also, \(f>t) = 
\4>i) = 0. These substitutions lead to tp = <f>. Plugging this into equations [14] 
and [D] we get 



JT 1 

2 <t> 2 (22) 



JL _ I _ 121 + 0(v 2 ) 
iv 2 24 ^ I 
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and 



p(n max ) = 1 - 0(v 2 ) (23) 



which agree with the results of the paper. If we examine equation [13| in this case 
we get 

p(n) = 1 - COS[( 2 2 " + 1) ^ ] (24) 
as the probability of measuring a target state after n iterations of Q. 



4 Application of the Formula for p(n) 

Next, we show how to apply our analytic formula for the probability of success 
after n iterations, P(n), to slightly speed up the quantum unstructured search 
algorithm. Although the speedup we obtain is not dramatic, it is worth making 
the point that it is possible at all as Christoph Zalka has proved, correctly, that 



Grover's algorithm is exactly optimal 25 . Many people have assumed, there 



fore, that it is impossible to beat Grover's algorithm. However, by combining 
techniques of quantum computing with those of classical computing we show 
that it is possible to do a little bit better than Grover's algorithm on average. 
The result we report was apparently discovered previously by Boyer et al. || It 
is shown here to persist for the case of fully generalized quantum search. 

We consider a punctuated quantum search algorithm that works as follows: 

Algorithm: Punctuated Quantum Search 

1. Run the quantum search algorithm for n iterations. 

2. Read the memory register. 

3. If the result is a target state halt, else reset the register to the starting 
superposition and return to step 1. 

The average time, T avg (n), it will take to find a target state if we stop the 
generalized quantum search algorithm after n iterations of Q is 

T avg {n) = YaLi (l-p(n)y~ p(n)in 

= ^TT (25) 

_ 2n 

— l-cos[(2n+l)0] ' 

We can find the optimal strategy, i.e., the best number of iterations to use before 
we attempt to measure the register, by minimizing the expected running time 
T a vg- To do this, we set the derivative of T avg to zero and solve for n — n opt . 

2-2 cos[(2n + l)<f>] - 4n</> sin[(2n + 1)0] _ 

(l-cos[(2n+l)0]) 2 

Typically n will be much larger than one so we can make the approximation 
(2n + 1)4> — 2n(p = x. By dropping the denominator we obtain 
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1 — cos x = x sin x 

2 sin 2 (f) = 2a; sin (f) cos (f) (27) 
x = tan (|) . 

which gives x — 2.3311 as the lowest positive solution. This solution corresponds 
to the minimum of the function. Hence the optimal value of n is 

n op t ^ — - — (28) 

<P 

This value of n corresponds to an average number of iterations of 

r l \ L38Q1 /0Q\ 

T a vg{n = nopt) - — 7 — ( 29 ) 

compared to 1,5 < J 08 iterations if we run Grover's algorithm until the probability 
is maximal. 

It is interesting to note that if we restrict the analysis some more to the 
case where \a) is given by equation ^, and there is only one target state then 
T aV g(ri opt ) ~ 0.6900a/F- This is faster than the lower bounds in Q, 0, §, 
and p5| ], but we are using a somewhat different model. They are looking at 
the minimum time it would take without measuring to find a solution with 
certainty up to errors from rounding n max to the nearest integer. The model 
we use allows for punctuated measurements and resets of the quantum search 
algorithm. Nevertheless, the punctuated quantum search algorithm is faster on 
average. Note that we have assumed that the time it takes to measure and reset 
the algorithm is negligible. This is reasonable as it only requires one function 
call. 

The punctuated quantum search algorithm has another advantage in that 
it should make it easier to eliminate decoherence. If we wait until we have the 
maximal probability of measuring a target state then we must maintain coher- 
ence for T steps as compared to only 1A ® 55 steps for the fastest measure 
and restart method. In fact if we are willing to settle for an average time equal 
to the time it takes to have maximal probability then coherence need only be 
maintained for 0,7854 steps at a time. 



5 fc-Parallel Quantum Search 



A way to speed up Grover's algorithm still further is to have a society of k 
computational agents all running Grover's algorithm independently at the same 
time. This is promising because the standard deviation 

a T = -^- J[l-p(n)}[l-p(n)+p(nr} (30) 
p{n) 

of the method we have already analyzed is fairly large and hence having multiple 
algorithms running may give a considerable speed up. 



9 



Suppose that we know that there are exactly r solutions amongst N can- 
didates. Given p(n,r,N), the probability of success for a single agent after n 
iterations, we can boost the success probability by using k agents acting in 
parallel. In particular, the probability that at least one agent, in a society of 
k independent agents, succeeds after each agent has undergone n iterations is 
given by 

Pk (n,r,N) = 1-(1 -p(n,r,N)) k (31) 
Thus the expected cost, Tavl, of performing fc-parallel quantum search is given 

by 

oo 

T iva = y 3 n Pk (n,r,N)(l- p fe (n,r, A0) i_1 = . - — = w 

aVS U Pk(n,r,N) l-(cos(I(l + 2n)arccos(l-^)) 2fe 

(32) 

As in equation ^ we can find the value of n that minimizes the expected cost. 
For <C 1, i.e., when there are very few solutions amongst the items searched 
over, arccos(l — ^) w 2^J^. Hence the average cost for fc-parallel quantum 
search is given by: 

T«(n,r, AT) « — ? — (33) 

l-(cos((l + 2n) v ^)) 2 ' £ 

gT^) (n r N) 

To find the mimimum, we find where — av9 g n ' — is equal to zero. This deriva- 
tive is given by: 

dT<tf g (n,r,N) = 1 - (cos((l + 2n)^)f k (l + 4fcn^tan((l + 2n)^)) 
dU (-l + (cos((l + 2n) v ^)) 2 ' c ) 2 

(34) 

Substituting x — (1 + 2n)y^ and realizing that n 3> 1 we obtain the 
following: 

dT^l(n,r,N) 1 - (cos(a;)) 2fe (1 + 2kxttm{x)) 

— (35) 



On 



(-i+tcos^y 



The variable x < 1 provided n < | (\J~^ ~ -0 • ^ e know that we can solve 
the problem with near certainty if we iterate Grover's algorithm to the maximum 
probability state in O(j^J^) iterations. Hence, for a large enough number of 
parallel search agents, k, there is a reasonable chance that the optimum number 
of iterations, n pti m um(r, N,k) at which the expected search cost is minimized, 
satisfies the criterion that x < 1. We therefore expand equation [55] as a series 
approximation in x about x — to order 0{x 2 ). Hence 

dT&l(n,r, N) _ -60 + 10(-1 + 3k)x 2 + 3(-l + 5fc 2 )z 4 
dn 60kx 2 
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As this equation is fourth order in x it can be solved analytically. Three 
of the roots are non-physical but one corresponds to an approximation to the 

true minimum of T^} g {n,r,N). Specifically, we find that dT <"'^' r ' N ' > = and 

Tavl(n,r,N) is minimized when x is given by 



1 5 - 15k + - 30k + 225F 

X ~ V -3 + 15fc 2 



(37) 



We note that x < 1 for all k > 2. Hence, the derivation of the optimum 
formula is self-consistent. Hence, as x = (1 + 2n)-J^, we obtain the formula for 
n op timai{r, N, k), the predicted optimal number of iterations to use for each of k 
quantum searches acting independently in parallel as 

. a r 7 \ W k- 15k + V5V-31 - 30fc + 225fc 2 
n optima i(r, N,k) » - ^ Zs + iEi? 

Hence the expected cost for optimal fc-parallel quantum search is given ex- 
plicitly by: 




/ 5-15fc+V5V-31-30fc+225fc 2: ' /~N _ ^ 

T£l(n,r,N) « — V lg±^ V (39) 

2 - 2cOS 2fc f y 5-15fc+vWgT=-30fc+225fc^ j 



6 Conclusions 



In this paper we have shown how to generalize the analysis of unstructured 
quantum search to incorporate the effects of an arbitrary starting superposition 
and an arbitrary unitary operator (or, equivalently, arbitrary averaging state) 
simultaneously. We have also shown that, rather than iterating the amplitude 
amplification operator until the maximum probability of success state it attained 
(i.e., for O(0.785398\/iV) iterations) it is better to stop after only 0(0.6900^) 
iterations (i.e., 88% of the maximum probability case). This strategy, is therefore 
approximately 12% faster than Grover's algorithm on average. 

Moreover, an ever better quantum search algorithm can be obtained by run- 
ning k independent quantum searches in parallel, stopping as soon as any of the 
quantum searches finds a solution. We find that the optimal fc-parallcl punctu- 
ated quantum search strategy is different from that of single agent punctuated 
quantum search strategy. In general, the higher the degree of classical parallelism 
the less (parallel) time is needed to perform the quantum computation. This intu- 
ition is captured analytically in equation (38), which gives the explicit connection 
between the number of amplitude amplification iterations as a function of the 
degree of parallelism k. This result is of practical utility to experimentalists. In 
particular, in any physical embodiment of a quantum search there will be some 
natural coherence time beyond which the computation becomes unreliable. Of 
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course, quantum error correction allows this time to be extended greatly, ar- 
guably indefinitively, if the individual error probability per gate operation can 
be made sufficiently small. While we believe this to be true, in practice it might 
be extraordinarily difficult to achieve. Instead, if we can predict the degree of 
parallelism needed so that the quantum search has a good chance of complet- 
ing within the natural coherence time of the physical system being used as the 
quantum computer, then the strategy of massive parallelism might provide a 
realistic alternative to relying solely on quantum error correction. Thus we see 
the classical parallelism as an adjunct to quantum error correction rather than 
a replacement for it. Equation (38) exposes precisely the space/time tradeoff 
between quantum coherent computing and classical parallelism, at least in the 
context of unstructured quantum search. 
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Fig. 1. Plot of the probability of success of Grover's algorithm after n iterations 
of amplitude amplification when there are r solutions amongst N = 64 possi- 
bilities. White regions correspond to probability 1, black regions correspond to 
probability 0. Notice the periodicity in the success probability as the number of 
iterations grows. 
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Fig. 2. Plot of the optimal number of iterations to use in /c-parallel quantum 
search as a function of the degree of parallelism k for r = 1 to r = 5 solutions 
(top to bottom in the figure) for the case of a database of size N — 2 20 . The 
solid curves are produced by numerical optimization. The points are given by 
our approximate formula for n ptimal ( r ; N, k) 



